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We develop a time dependent random matrix theory describing the influence of a time-dependent 
perturbation on mesoscopic conductance fluctuations in open quantum dots. The efi'ect of external 
field is taken into account to all orders of perturbation theory, and our results are applicable to both 
weak and strong fields. We obtain temperature and magnetic field dependences of conductance 
fluctuations. The amplitude of conductance fluctuations is determined by electron temperature in 
the leads rather than by the width of electron distribution function in the dot. The asymmetry 
of conductance with respect to inversion of applied magnetic fleld is the main feature allowing to 
distinguish the effect of direct suppression of quantum interference from the simple heating if the 
frequency of external radiation is larger than the temperature of the leads hu> ^ T. 

PACS numbers; 73.23.Ad, 72.15.Rn, 72.70.+m 



Introduction - Transport coefficients of disordered 
and chaotic electron systems fluctuate from sample to 
sampleErQ. These fluctuations are commonly called 
mesoscopic fluctuations. Mesoscopic fluctuations of con- 
ductance of non-interacting systems are universal. The 
universality means that the variance of the conductance 
((5(7^) is of the order Gq, where Go = e'^/nh is the quan- 
tum of conductance and is weakly dependent of the sam- 
ple geometry, seeD'El. 

The fluctuations of transport properties of electron sys- 
tems is a quantum mechanical phenomenon based on 
the interference of quantum states. As any other in- 
terference phenomena, conductance fluctuations are very 
sensitive to inelastic processes, commonly referred to as 
dephasingp. The dephasing processes in open quantum 
dots -were considered on the purely phcnomcnological 
basiaj. First microscopic consideration of the microwave 
radiation on the weak localization in quantum dots was 
performed in Ref.tJ. In this reference, the concept of 
the time dependent random matrix theory (TRMT) was 
used. 

The purppsj3 of the present paper is to extend the re- 
sults of Ref.ll3 to describe the effect of the external mi- 
crowave radiation on the mesoscopic fluctuations of the 
conductance. The ultimate goal is to identify observable 



(1) 



features which allow one to distinguish the effect o; 
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the 



external field to the dot itself from a simple heatingZl. 

However, there is a significant difference in the calcula- 
tion of the mesoscopicj-conductance fluctuations and the 
averaged conductancell3. Because, the dot is subjected 
to the external classical radiation which produces non- 
equilibrium in the dot, the d.c. -current Iq through the dot 
is finite (though randomly changing from one configura- 
tion to another) even if the d. c.-voltage V = Vi — V^ across 
the dot is zero, see Fig. 1. Kiis current /q is due either to 
the photovoltaic effect (seetJ and references-Lherein) or to 
rectification of a. c. -bias across the dot, seell3. We are in- 
terested in the linear response to the applied d. c. -voltage 
across the dot: 



In principal, the linear in V contribution to the current 
comes comes from two sources; (i) the non-equilibrium 
of the distribution functions in the leads, (ii) change in 
the photovoltaic current, correspondent to a different re- 
alization of the dot due to the finite bias. Nevertheless, 
we will show that due to the electrOpueutrality condition 
the non-equilibrium current prevailslia. 

Closing the introductory part, we note that the present 
paper has a certain pwerlap with the recent preprint by 
Wang and Kravtsovcil, where the conductance fluctua- 
tions were calculated for open quantum dots subjected 
to a periodic ac pumping. Our treatment is different 
in several aspects. Firstly, our results are applicable for 
the frequencies of the external radiation uj smaller than 



the T: 
of Ref. 
like Ref. 



iless energy of the dot, Et, whereas treatment 
valid in the opposite regime. Secondly, un- 
Lil, we will restrict ourselves to the case of the 
monochromatic radiation acting on the dot. Finally, we 
will highlight the role of the electro-neutrality require- 
ment in a separability of the photovoltaic effect and the 
mesoscopic conductance fluctuations, which was not done 
in Ref.y. 

Model — We apply the random matrix theory (RMT) 
to study the conductance of open quantum dots, see 
Ref.Q. All corrections to the RMT are governed by a 
small parameter A'ch/ffdot, where gdnti = Et/Si and 6i 
is the mean level spacing, see Ref.El for the detailed 
discussion. We consider the conductance fluctuations 
of quantum dots with a large number of open channels 
A^ch 3> 1. In this approximation, we neglect the effects 
of the electron-electron inteeaetion on the conductance 
which are as small as I/A^^j^EjO, while the conductance 
fluctuations are proportional to 1/A^ch- The same condi- 
tion also iflJlows us to use a conventional diagrammatic 
techniquetj to take the ensemble average. External mi- 
crowave radiation is modeled as time dependent random 
part of the Hamiltonian of the dot. 
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FIG. 1. An open quantum dot is connected to two leads 
with applied voltages Vi,r- The measured current through the 
dot has an offset Jo at zero bias and the linear response to 
small applied voltage V — Vi — Vr- 
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The Hamiltonian of the system i 



(2) 



where Hb is the Hamiltonian of the electrons in the dot, 
which is determined by the AI x M matrix Hnm- 



M 



(3) 



4>n corresponds to the states of the dot and the thermo- 
dynamic limit M — ^ 00 is assumed, Ec is the charging 
energy, and n = J2m=i the last term in Eq. (||) 

is the largest contribution to the interaction effects in 
quantum dot, see RefJlJ for the discussion of the status 
of this approximation. Matrix HY){t) is given by 



Huit) =H + Vip{t). 



(4) 



Here the time independent part of the Hamiltonian H is 
a random realization of a ill x M matrix, which obeys 
the correlation function 



(i/„„(ci>i)i7:,„,(ci>2))=Af 



1 - ^) Snn'S., 
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where 5i is the mean level spacing of the dot and pa- 
rameters iVrfTc describe the effect of the magnetic field 
on the dotEj. These parameters can be estimated as 
iVd,c =i 5dot (*i T ^"2)^ /*o where $1,2 is the magnetic 
flux through the dot and $0 = hc/e is the flux quantum. 
The time dependent perturbation is described by sym- 
metric M X M matrices Vnm and function ip{t) of time 
t. We assume that the perturbation is harmonic with 
single frequency w, ^p{t) = coswi, even though most of 
the consideration [up to Eq. (p7|)] is valid for an arbitrary 



function. The effect of the perturbation-pn the system is 
totally determined by two parameters,E3 



(6) 



Parameter Z has a meaning of the average velocity of 
the energy levels of the dot under the external pertur- 
bation and can be omitted from our consideration due 
to screening [see belcwd . The parameter Co characterizes 
its typical deviation|13. Since all the physical responses 
of the system are characterized by the same parameters, 
the value of Co can be eliminated by an independent 
measurement. 

The electron spectrum in the leads near Fermi surface 
can be linearized: 



(7) 



where ipaik) denotes different electron states in the leads, 
k labels the continuum of momentum states in each chan- 
nel a, hvp = l/27rt^ is the Fermi velocity and is the 
density of states per channel at the Fermi surface. We 
put ?i = 1 in all intermediate formulas below. 
The coupling between the dot and the leads is 



(8) 



For the reflectionless point contacts—the coupling con- 
stants, Wnai in Eq-(||) arc given byotj: 




if n a < iVch, (9) 
otherwise. 



For open dots with a large number of open channels 
-^ch 3> 1 the interaction term can be treated within mean 
field approximation, so that the Hamiltonian takes 
the form 

M 

n'S'{t)= J2 ^i[HDnm{t)+eVdit)Snm]4'm, (10a) 
n,m— 1 

eVd = 2E,{n)g, (10b) 

where {n)g stands for the quantum mechanical (but not 
ensemble) of the number of elec tron s in the dot. Cor- 
rectiops to mean-field treatment ( |lOa[) were calculated in 
Ref.tj and shown to be small as 1/ NX- 

In the mean field approximation (|lOa| ), one can intro- 
duce one-particle S - matrix Sa0{ti,t2) as 

SMt, t') = 5o,p5{t ~ t') ~ 27nvWl,Gfi{t, t')W^p, (11) 

and the Green functions &^'^\t, t') are the solutions of: 

d 



= S{t-t'), 



(12) 
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where the matrices Hd and W are defined by Eqs. 
and (^). 

The d.c.-current through the dot is given by, se 



(s^{t2MS{t, h) - k5{h - h)) } (13) 



where Tp is the period of the external perturbation, f[t) 
is related to the Fourier transform of the electron distri- 
bution function in the ath channel as 



iTe 



sinh nTt 



(14) 



and 



A, 



Q/3 — Sap < 



iVch' 



for 1 < a < 



for M + 1 < a < iVch- 



(15) 



The spin degeneracy is taken into account in Eq. (13) 



We assume that the degeneracy is not lifted by magnetic 
field. 

To complete the theory one needs an equation for the 
averaged number of particles {n)q, see Eq. ( |10b| ). It is 
found from the continuity relation as 



d{n{t))g 
dt 



dtidt2^T: 



[f{tl - t2) 



(16) 



(s\t2,t)S{tM)-5{t2-h))] 



Eauations (|T|) — (|l|) are similar to those used in 
Ref.cj for studying the frequency dependence of the con- 
ductance of mesoscopic system. 

Ensemble averaging — Our goal now is to perform cal- 
culations of the conductance correlation function 



7^($l, $2) - (5($i).9($2)) - (g(1>i))(ff(1>2)), (17) 



using the model outlined above. 




• -f(t-t') 
t' 

FIG. 2. Two diagrams, which contribute to the conduc- 
tance correlation function R. 

We use the leading approximation in small parameter 
1/^ch- The fluctuations of the conductance are smaller 
than its average value and we can use instead of sample 
specific Eq. (Hq) its ensemble averaged counterpart: 

d{n{t))q _ N^iA 



dt 

eNi 



2ti 



(18) 



eNr 



+— (M - - Zip{t)) + —L(y,^VA- Zip{t)) 



Equation (|18|) is nothing but a discrete form of the diffu- 
sion equation for the bulk system and the last two terms 
correspond to the divergence of the drift current. Sub- 
stituting Eq. ( |lOh| ) into Eq. ([l8|), solving the resulting 
differential equation, we find 

AeE^ NxVi + N,V, 



eVS) + Zip{t) = 



4£;c + (Si TVch 

(5i + {2^/N,^)dt 

Si 



(19) 



Z- 



iEc + (27r/iVch)at 

We notice from Eq. (^^ that the characteristic energy 
scale governing charge dynamics is EcNch/^n. Usually, 
this scale is of the order of the Thouless energy, Et- Be- 
cause all the random matrix theory is capable to describe 
the energy scale only smaller than Et, we can consider 
only uj <C Et — EcNch/^ir. Moreover, for the small 
quantum dot Ec^ 5i, so that Eq. (19) gives 

mvi + N.v. ^^^^ 



ch 



and the time dependent perturbation (|^) can be consid- 
ered as traceless, Z = 0. This constant in time compo- 
nent of the bias of the dot can be removed from Eq. (02) 
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for the Green function by the following gauge transfor- 
mation: 



G(t,t')^G (t,t')e 

Vd=0 



'\„-ieVd{t-t') 



(21) 



Substituting Eq. ( |2l] ) into Eq. ( [T^ ) and expanding up 
to the first power m V = Vi — K-, we find 



5ci — Go 



NiN, 
Neb 



is the classical conductance of the dot, Go — /ttTi is the 
quantum conductance and F{x) is the Fourier transform 
of the derivative of electron distribution function: 



dl 



gy - 9ci + Go jdt j dtidt2F{T{ti (22) 

-oo 

X Tr|5(i,ii)A5^(t2,t)A|. 



Here 



F{x) 



sinh; 



(23) 



The correlation function of the conductance ([1^) is 
given by the diagrams shown in Fig. 2 and can be found 
from the following analytical expression 



i? = 2 



3cVf [dtdf 



47r2 / T2 



OO OO 

jdrF^iTr) j dO 



V h 0, ,t-t']vi h 0, h -, t - t 

2'2 2' /l2'22' 











r/2 



(24) 



The diffuson and the Cooperon are defined by the fol- 
lowing equations: 

C(Ti,T2,i) = e(Ti-r2)exp(^-i^ \c{T,t)dT^ , (25a) 



I?(ti,i2,r) = e(ii-i2)exp 



where 



Kd{T,t)dt], (25b) 



Kd,c = 7d,c + Co {ifiit + r/2) - ip{t - r/2))' (26a) 



7d,c = TT i^ch + N^.c) , 

ZTT 



(26b) 



and parameters A^d.c describe the effect of the magnetic 
field, see Eq. d). 
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Derivation of Eq. ( P6b| ) deserves a little bit of addi- 
tional discussion. One notices, that the diagrams of Fig. 3 
does not contain any piece corresponding to the classical 
distribution function in the dot. We can trace it into the 
expression for conductance (^) which contains traceless 
vertices A, which can not be dressed by the dashed line. 
On the other hand, any vertex with finite trace corre- 
sponds the modified distribution function oLelectrons in 
the dot and represents the effect of heating EiJ. Since the 
distribution function is not dressed in the expression for 
conductance fluctuations, we conclude, that the effect of 
heating is not relevant for the conductance fluctuations 
and the temperature dependence of the conductance fluc- 
tuations is uniquely determined by the electron temper- 
ature in the lead£3. That meaafi that contrary to the 
common believes, see e.g. Ref.Q the amplitude of the 
mesoscopic fluctuations can not be used for the study of 
the distribution function of electrons in the dot. From the 
theoretical side, it is important to emphasize, that the ap- 
pearance of the traceless vertices is determined solely by 
the electro- neutrality condition (|^, any other choice of 
the dot bias would lead to the change in the photovoltaic 
current til. 

Limiting cases — Below we consider the limit of high 
{huj 3> C) and low {huj <^ jd,c) frequencies. For the high 
frequencies, w S> C, we obtain 



(c) 





rco 
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V7d 


7d/ 



^Ti"' V 7d V 7d ' 7d / 7c V 7c ' 7c ' 7c / 



(27) 



(d) 

FIG. 3. The elements of the diagrammatic technique. where dimensionless Q— functions are given by 
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{x,y)= / drF^iyr) 



(28a) 



27r 



exp (-T (l + 2xs\n^ C,/2)) dC, 
l + 2xsm^ C/2 2^' 



Qc{x,y,z) = / drF'^iyr) 



1-K 



exp (-T (1 + 2xsiii^ C/2)) 



(28b) 



1 + a;(sin^ C/2 + sin^(C/2 + zt/2)) 27r 



Let us now consider the dependence of the functions 
Qd,c and Qc{x,y, z) from Eq. ( 28b ) on temperature y. 
For the limit of high temperature, y ^ 1, we obtain 



zix,y,z) « Qd{x,y) 



1 



3y VI + 2x 



(29) 



The equahty between functions Qc and (5dr*^ieans the 
magnetic fi eld symmetry of the conductanceEil. Indeed, 
using Eqs.( ^6a|) , ([2^ ) and ( ^ ) we observe, that 



i?($i,$2) - i?($i,-$2) 



(30) 



However, in low temperature limit y <C 1, wc obtain for 
X > 1 



Qd{x,0) 



1 + x 



1 



(l + 2x)3/2 
Qc{x,0,z) w 



(31a) 
(31b) 
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FIG. 4. Functions Qc{x,y,z) and Qd{x,y) computed for 
X = 1 and z = 10. As temperature y — T/7 increases, func- 
tion Qc{x,y,z) approaches frequency independent function 
Qd{x,y). 



Comparison of Eqs.(26a), ( |27| ) and (31a) reveals an 
important fingerprint of the dephasing by the external 
radiation — violation of the Onsager relation 



i?($l,$2) 



Co' 



(32) 



where 7 = 7d('i'i,^'2) = 7c(^i, ~*i'2), provided that 
7c(<i>i,$2) = 7d(<i>i, -<i>2) > 7- This breakdown of the 
Onsager relation is a simple manifestation of the lifting 
of the time reversal symmetry in the system with time 
dependent Hamiltonian. 

In the limit of low frequency huj <^ 7d,cj the contribu- 
tion from the Cooperon and diffuson parts are described 
by the same function, so that the conductance correlation 
function can be represented in the form 



R = 



47r2 



7d V7d Id J 7c V7c 7c 



so the Onsager relation ( pO[ ) holds. Here, 

r^'' d^dC 
lo 



Q{x,y) 



47r2 



(33) 



(34) 







X 



exp(-(H-4xsin2 ^/2sin2 C/2)t) 
1 -t-4xsin2^/2sin2 C/2 



dr. 



This expression in the limit of high temperature T ^ 7d., 
has an asymptotic behavior 



Q[x,y)^—K{-^x). 
32/ 



(35) 



At zero temperature Q(x, y) is given by the expression 

(36) 



l i;(-4.T) + (l+4:r)j^(-4a:) 



where K{x) and E{x) are the elliptic integrals of the first 
and second kind respectively 



K{x) 



7r/2 



dip 



'0 \/\ ~ X SVC? ip 

E{x) = / \J \ — xsiv? ip dip. 
Jo 

We conclude that the conductance fluctuations are 
suppressed by external radiation even in the limit of the 
low frequency, see Eq. (^3|). Indeed, during one period of 
time, the system goes along a closed loop in the param- 
eter space and the contribution to the d.c. -conductance 
is effectively determined by the equilibrium conductance, 
correspondent to each point of the loop. The equilibrium 
conductance fluctuates along this loop. Thus, the ob- 
served d.c. -conductance is already partially averaged over 
some realizations of the quantum dot and its fluctuations 
decrease. The perturbation strength is related to the 
length of the contour in the parameter space and effec- 
tively determines how many different dot's configurations 
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contribute to the d.c- conductance. Consequently, the 
stronger perturbation, over the larger number of the real- 
izations the d.c. -conductance is averaged and the smaller 
fluctuations of the d.c. -conductance. 

This should be contrastecLwith the suppression of the 
averaged magnetoresistancaS'tJ. There, the stationary 
field does not do anything because the result is already 
ensemble averaged. In order to suppress the average 
quantum correction, the field should have change on the 
time scale of the order of l/7esc, where 7osc = SiNch/^ir 
is the escape rate from the dot. That is why the effect of 
the low-frequency radiation on conductance fluctuations 
and weak localization corrections are significantly differ- 
ent. At high frequency huj ^ 7esc the d.c. -conductance 
no longer can be represented in terms of the stationary 
conductance and the suppression of both the conductance 
fiuctuations and the weak localization correction to the 
conductivity can be interpreted as dephasing. 

Comparison with experiment — Our results still con- 
tains an unknown parameter Cq characterizing the 
strength of the perturbation. There is a way, however, 
to present the results in a form not depending on this 
parameter, thus eliminating a need for additional fitting. 
Following Ref.Q, we represent the parametric dependence 
of the weak localizatipa correction (5g„i versus var g, 
where Sg^ii is given bytEl: 



0.3 



gwi 



-P 



Co huj 

Tcsc 7csc 



P{x,z) 



nh {Ni + iVr)2 



sin 



(37) 
(38) 







The conductance variance is determined from Eq. ( |24| ) 
for broken time-reversal symmetry 7c — > oo. Figure 5 
shows the parametric dependence for various values of 
the parameters Co and lo and T = 107osc- 

We observe that the shape of the curves depends on the 
frequency of external radiation. Particularly, in the limit 
of low frequency huj <^ j^sp the weak localization correc- 
tion is not changed by the radiation, while the conduc- 
tance fluctuations may be significantly suppressed. At 
high frequency, huj 3> Cq, 7csc the curves become non- 
sensitive to the radiation frequency. 

The authors of Ref.i found that the radiation applied 
to their device produced curves in varg vs. dg^/i plane 
identical to the curve produced by increasing tempera- 
ture of the device for a wide range of frequencies. This 
observation apparently demonstrates that the radiation 
produces the heating of electronSp-and the effect of de- 
phasing without heating, see RefJ23, is not observed in 
experimentaj. It was also suggested that the main mech- 
anism is the increase of the temperature in the dot due 
to the Joule heat by induced ac source-drain bias. 
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FIG. 5. Weak localization correction Sg„i versus con- 
ductance fluctuations var g of an open quantum dot with 
A^i = A'^r 3> 1 for three values of frequency u: Tiuj = 0.57csc 
(o); hio — 57csc (v); hto = 507esc (A). The temperature for 
all lines was taken T = 107esc. The amplitude of the field Co 
varies from 10~^7csc to 10^7csc. Go ~ /-kTi is the quantum 
conductance. 

Although the present data of Ref.0 support the heat- 
ing mechanism of suppression of the weak localization 
correction to the conductance and the conductance fluc- 
tuations, we believe that more detailed analysis has to 
be done. According to our theory, see Eq. ( p2| ) and the 
paragraph below Eq. (26b), (i) mesoscopic fluctuations 
are sensitive only to the temperature in the leads, and 
therefore, the notion of heating of electrons in the dot 
responsible for 1 /T dependence of the mesoscopic fluctu- 
ations is not relevant: if there is a heating, it manifests 
itself only through the temperature of the leads; (ii) high 
frequency curves of our theoretical Fig. 5 quantitatively 
agrees with data on Fig. 3 of Rcf.El, for frequencies / = 1, 
10 and 25 GHz. An exception is the lowest frequency 
curve (/ = 100 MHz) represented in this plot, for the dot 
with Si — 2.4/iey, iVch = 2 corresponds to hf j^esc ~ 0.5, 
so according to our Fig. 5 it should have observable devi- 
ations from the high frequency curves, which is not seen. 
However, taking into account uncertainty in determina- 
tion of the level spacing 6i from the geometrical area of 
the dot, this does not unambiguously rule out the mi- 
crowave dephasing mechanism. 

We believe that the "smoking gun" evidence for the 
mechanism considered in the present paper is the vi- 
olation of the Onsager relation ( ^ ) in high frequency 
regime, Ulj ^ 7c, d, T . The dependence of this violation 
on the amplitude of the fleld Cq is the main prediction 
of the theory. 

Summary — We constructed the time dependent ran- 
dom matrix theory to describe the effect of the non- 
equilibrium external radiation on conductance fluctua- 
tions of an open quantum dot. The main experimen- 
tal feature to reveal such a mechanism is the breakdown 
of the magnetic symmetry of the conductance by high- 
frequency radiation, huj ^ T, see Eq. (p2[). 
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